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Abstract: In this paper,we investigate the mathematical modeling for the cosmological attractors 
propagated in mimetic gravity upon which an interacting dark energy-dark matter is supposed to 
be existed. The average value of the interaction of these percentages, namely T, say, may be used 
to investigate generally the modeling of an attractor; the actual value could only be determined by 
data in any particular case. We have seen, for example, that it was led to investigate the subject of 
initially invariant submanifolds. 
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Introduction In many passages of our studies on large 
scale structure of the whole Universe, observational data 
had distinctly stated that the current epoch of Universe 
which it dominated here give an explanation of the accel¬ 
eration expansion of the Universe. Namely the data from 
supernovae SNe la [Tj, cosmic microwave background 
radiations via WMAP [2], galaxy redshift surveys via 
SDSS [3] and galactic X-ray gj. But only the means 
of observing and connecting them together with theory; 
that the meaning of all phenomena, and the reason of 
their peculiar connexions, was a philosophical problem 
which required to be attacked from a different point of 
view; and that the significance especially which lay in 
the phenomena of matter contents and exotic form of 
it would only unfold itself if by an exhaustive survey of 
the entire fully theoretical description, individually, con¬ 
sistently, and physically, we gain the necessary data for 
deciding what meaning attaches to the existence of this 
current acceleration expansion, or dark energy (DE) and 
dark matter (DM) in the contents of the universe. DE 
is considered as an exotic matter with negative effective 
pressure to energy density ratio in a very recent epoch 
z ~ 0.7. The existence of such type of energy brings 
two important unsolved problems: The fine tuning and 
cosmic coincidence. Perhaps the best candidates for DE 
are the quintessence scalar field or a phantom energy field 
[5] . The geometrical modifications of Einstein gravity are 
also interesting and have been studied widely in litera¬ 
tures 0 -na. In the case, however, of many observatories 
m -(19 , especially as regards the older records, many 
data for interacting DE-DM exist; further, the decay of 
DE to the DM and unparticle is at best only an approx¬ 
imation [20].|21j. the success attending which probably 
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varies considerably at different stations [22 ■ When the 
continuity equations are satisfied for different types of 
matter(energy) contents with interaction terms, a phase 
space description is found to be useful [23]-[25]. In this 
approach we present, a suitable multidimensional rep¬ 
resentation of a first order dynamical system in which 
each dimension corresponds to one density of the system. 
Thus, a point of phase space corresponds to a specific 
state of the system, and a path represents the evolution 
of the density functions of the system through different 
states. Our aim in this letter is to perform phase space 
description for an interacting DE-DM system in a type 
of modified gravity, called as mimetic gravity (MG). 

Basic equations Another opportunity to solve DE 
problem will present a type of disformal invariant 
theory, where we might be identified with a specific 
reparametrization of the metric = (g a ^daCpdpejfjg^v 
[26]: 

S = ~\J d*xy/-g(~g^) +Cmh (1) 

Here we suppose that g^ represents an auxiliary metric 
and 4> denotes an auxiliary (non ghost)scalar field. Such 
type of modified gravity recently has been studied from 
different points of view [27]-[35]. Derivation of the field 
equations are straightforward. If we vary the action ([I]) 
with respect to the metric g and the scalar field we 
obtain the following set of field equations: 

(G^ - k 2 T^) - (G - K 2 T)g> lo ‘g vp d ol 4)dp(j) = 0, (2) 

A=9 k [^(G - n 2 T)g* x d x ct>\ = V K [(G - T)3 K </>] = 0, 

(3) 

The acceleration expansion of Universe in standard cos¬ 
mology may be attributed to a homogenous and isotropic 
FLRW metric in the following form: 

ds 2 = dt 2 — a(t) 2 dx 2 . 


(4) 






2 


The set of FLRW Eqs. 
metric Q: 


can be obtained using (2|3) for 


3 H 2 = K 2 p e ff = K 2 ( 3 Ptot - ( 5 ) 

R 

2 H = -K 2 p e ff == k 2 ( 4 Ptot - - 2 ) ( 6 ) 


here p tot = E f =1 Pi = Ef =1 w i p i . 

We assume the matter contents are composed by three 
component fluid containing matter p m , dark energy pd 
and radiation p r . We suppose that these components are 
interacting with an unknown set of interction fuctions IV 
The corresponding continuity equations for three compo¬ 
nents can be written in the following system: 


Pd + 3H(pd+Pd) — Ti, 

p m + 3Hp m = T 2 , (7) 

Pr + 3H(p r + Pr) = r 3 , 


To preserve the total continuty, the interaction functions 
must satisfy collectively such that ^+^ + ^ = 0. To 
define an autonomous dynamical system , it is needed to 
define dimensionless density parameters via the following 
set of new variables: 


_ K 2 p d _ H 2 p m ~ _ K 2 p r 

“3 H 2 ' V ~ 3H 2 ’ ^ “ 3 H 2 ' 


Using the FLRW Eqs we can calculate the following im¬ 
portant quantity: 


for any practical purpose to which the above systems 
would be applicable in the ordinary analysing of a local 
stability point. But this was based upon the assumption 
of a density-TL relation between interaction terms, the 
interaction function of which depended on their rela¬ 
tive densities pi as well as on Hubble parameter H. This 
being the case, we are at liberty to make the assumption 
that the interaction function of each term in ( To][l2 ) (un¬ 
der specified conditions) is known, without thereby in¬ 
troducing any risk of self-contradiction in mathematical 
calculations. This assumption has the great advantage, 
that the last interaction terms in ( Tofl2 ) now appears as 
the autonomous term of the components ( x,y,z ). Four 
models are specified for our dynamical system which we 
will study the local stability points using them [231 . 

Interacting model - I The first model has the following 
interaction terms: 


T 1 = -6bHp d , r 2 = F 3 = 3bHp d , 


(13) 


The model to be considered in a scenario in which dark 
energy will decay into matter and radiation. The param¬ 
eter b > 0 represents decay rate. 

Using (13), the system ( Top2 | we obtain: 

^ = -3x(w d - 2 (w d x + w r z)^j - 6 bx, (14) 

= 6 y{w d x + w r z ) + 3 bx, 

= —3 z[w r — 2 (w d x + w r z)^j + 3 bx, 


H 

Jr 2 


-3(s^ =1 w 0 a; 0 + i). 


(9) 


here x a = ( x,y,z ). We change time coordinate from t 
to the N = In a. Using we rewrite the continuity 
equations |7j) in following autonomous system in the di¬ 
mensionless variables: 

= -3x(w d - 2{w d x + w r z)^ + ^j-, (10) 

^ = 6y(w d x + w r z) + 1 ^^-, ( 11 ) 

^ = -3z(w r - 2(w d x + w r z)') + (12) 


we shall assume w m = 0, w r = ^ and — 1 < w d < 

The coupling functions Tj, i = 1,2,3 are considered 
as phenomenological general functions of the energy den¬ 
sities pi and the Hubble parameter H i.e. Ti(Hpi). If 
we forego this assumption, the form of interaction term 
is Tj ~ E bHp or ~ E H~ 1 p a p n . In strictness the 
interaction term T must be supposed to act upon the co¬ 
ordinates ([8]) in its actual condition, whereas in ( Top2 |, 
previously cited, the system is supposed to be absolutely 
autonomous. 

Analysis of stability in phase space In this phase space 
approach it is assumed that the system ( T(i|[T 2 1 is au¬ 
tonomous, and this assumption is sufficiently accurate 


The stationary or critical points are the solutions by 
equating the left hand sides of (30 


31) to zero. That 
assumption would create the possibility of the existence 
of stable point, that is, of any critical point of the sys¬ 
tem. The possible solutions and the eigenvaluesA of the 
linearized system ^ = Ay near these points (i.e. when 
X = X C + Y,\\Y\\ 

<C ||V C ||), as the solutions for the vector equation AV = 
AU, are classified as the following: 


• Ai: X c = (0,0,0), A i = {0, —3/2uv — 3/2u>d — 3 6 — 
3/2\ — 2b — w d + w r \, —3/2w r — 3/2wd — 36 + 3/2| — 
2 b - w d + w r |} 


• B V .X C = (0,0,1/2), A, : = {1,1 — 3/210^ — 36 — 
3/2 \w d + 2 b\, 1 — 3/2 w d — 3 b + 3/2\w d + 2 b\}. 


• C7 i: 


x c 


He 

Zc 


—w r Wd + Wd 2 + 4 bw d — 2 bw r + 4 b 2 
2(—w r w d + w d 2 + 2 bw d — bw r ) 
b (2 b — w r + w d ) 

2{—w r Wd + w d 2 + 2 bwd — bw r ) 
_ b ( w d + 2 6 ) _ 

2{—w r Wd + w d 2 + 2 bw d — bw r ) 


(15) 

(16) 
(17) 


A i = {3wd + 6b, 3wd + 6b, 6 6+3 w d — 1}- The follow¬ 
ing theorem helped us to classify the critical points: 
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FIG. 1: Model I: Solutions x(N),y(N), z(N) for Wd 

—0.7, b = 0.5. Solid x(N), dot y{N) and dash z(N) 



FIG. 3: Model I: Effective EoS w e ff = f or _ 

J J x ~ry~r z 

—0.7, b = 0.5. 

increased. Field equations can be written in the following 
forms: 



dx 

dN 

dy 

dN 

dz 

dN 


—3x(wd — 2(wdX + w r z)j — 3bx, (19) 
6 y(wdX + w r z) + 3 b(x — y), 

—3 z(w r — 2 (wdX + w r z ) J + 3 by, 


When we solved the Eqs. ^ = 0, we obtained all critical 
points, eigenvalues of the linearized system: 


FIG. 2: Model I: Phase space for Wd = —0.7,6 = 0.5. It 
shows an attractor behavior. 


• A 2 : {x = 0,y = 0,z = 0j, Ai = {-1, -36, -3(w d + 
b). It may be considered in the system as a stable 
point. 


Theorem : The point X c is asymptotically stable 
for the system ( Io][l2 ) if the real part of every eigen¬ 
value is negative. It is unstable if any eigenvalue 
has a positive real part. 

We conclude that the points A\,B\ are unstable, 
while the point C\ can be a stable (unstable) de¬ 
generate node. 


• B 2 : {x = 0, y = 0, z = 1/2},Ai = {1,1-36,1- 
3 (wd + 6)}. This point is unstable. 

• C 2 : {x = 0 ,y = —(3/2)6 + 1/2,2 = (3/2)6}, 
A j = {36, 36 — 1, —3^}. This critical point is also 
unstable. 

• D 2 : 


The density functions functions (x, y, z ), phase portrait 
and the effective EoS are shown numerically in the follow¬ 
ing Fig.( |l|2|3p , in which the parameters are taken as the 
Wd = —0.7, 6 = 0.5; the solutions for density functions 
are those determined by numerical algorithms, and the 
effective EoS w e ff f ollow by calculation in the manner 
explained for (30]|3l). 

Interacting model - II The second model of interaction 
is given by the following terms: 


w d - 3 - 3bw d 

6wj + 26 — 2 Wd ’ 

—6 + 3bwd + 36 2 
6 w% + 26 — 2 Wd ’ 
36 2 

6wj + 26 — 2 Wd 


( 20 ) 


r x = —3bHp d , T 2 = 3bH(p d -pm), T 3 =3bHp m . (18) 
The DE was a loss anyway, but the radiation would be 


Ai = 3w d , A 2 = 3(w d + 6), A 3 = —1 + 3(wd + 6). (21) 

This critical point can be considered as unstable if 
w d < -6. 
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FIG. 4: Model II: Solutions x(N),y(N), z(N) for w d 
—0.7, b = 0.5. Solid x(N), dot y(N) and dash z(N) 


FIG. 6: Model II: Effective EoS w e ff = Wd * + ™ rZ for Wd 

J J x~ry-\-z 

— 0.7, b = 0.5. 



FIG. 5: Model II: Phase space for Wd = —0.7,6 = 0.5. It 
shows an attractor behavior. 


Numerical solution, is obtained by imposing suitable 
IC and parameters into (191 . In Fig.(??) we see an 
increasing form for radiation while the DE is decreased 
as well as DM. 


An cosmological latetime attractor is detected in Fig. 

(??). 

Effective EoS is computed numerically. As we observe 
in Fig. (??), the EoS w e // > —1- So, the model can be 
considered as a cosmological ciable model for interacting 
DE-DM. 


Interacting Model - III A third model which we’ll con¬ 


sider is a scenario in which DE is decayed but DM and 
radiation are increasing: 


Fr = -6 bn 2 H 1 p d p r , T 2 = F 3 = 3 bn 2 H 1 p d p r ■ (22) 


The appropriate form of the dynamical system is taken 
as follows: 


dx 

dN 

dy 

dN 

dz 

dN 


—3 xyui d — 2{w d x + w r z)j — 18 bxz, (23) 
6 y(w d x + w r z) + 9 bxz, 

—3 z(w r — 2(w d x + w r z )) + %xz, 


The stationary points and the corresponding critical 
points for (23) are given as the following: 


• A 3 :{0,0, 0}, A,; = {0,-1,— 3d}. This point is un¬ 
stable. 


• B 3 : {x = 0, j/ = 0, z = 1/2 },Xi = {-3w d + 1 - 
96,1,1}. This point is also unstable. 

• C 3 : {x = 1/2 ,y = 0 ,z = 0},Ai = {-1 + 3w d + 
(9/2)b,3w d ,3w d }. This critical point can be stable 
if —1 + 3w d + (9/2)6 < 0. 

• D 3 : 

* =_ ~ 3wd + 1 ~ 96 (24) 

8 1 62 + 546^-96’ 1 ’ 

— 18u^ — 81bw d — 2 + 12w d — 816 2 + 276 

V = 1626 2 + 1086u>d - 186 ’ 

_ —2 w d + 6 w 2 d + 9bw d 

546 2 — 66 + 36bw d 
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FIG. 7: Model III: Solutions x(N),y(N), z(N) for w d = FIG. 8: Model III: Phase space for w d = -0.7,6 = 0.5. It 
—0.7,6 = 0.5. Solid x(N), dot y{N) and dash z(N) shows an attractor behavior. 


and eigenvalues: 

Ai = 
A2 = 
A3 = 

here 


3 w d 


96 + 6 Wd — 1 ’ 

9 bwd + v 7 3 A 
66 (9 b + 6 Wd — 1) 
9 bwd — \/3A 
66 (9 b + 6 w d — 1) 


= -1377 bw d + 4860 b 2 w d + 2592 w d 2 b 


A^ 

bw d 

+432 w d 3 + 96w d — 360 w d 2 — 129 6 6 2 + 180 6 — 
There is (un)stablity for this critical point. 



In Fig. ([ 7 ]) a numerical solution is developed for density 
function(s). All types of the densities are decaying. 

An attractor which was started from an initial point 
is observed in Fig. (??). An affective EoS with range 
w e ff < —0.5 is drawn in Fig. |9]). 

Interacting Model - IV 

An alternative form of interaction is defined by the 
following: 


Ti = -36 k 2 H 3 p d p r , 

T 2 = 36k H (PdPr PmPr^i 

r 3 = 3bn 2 H~ 1 p m p r . (29) 

The corresponding dynamical system is written in the 


FIG. 9: Model III: Effective EoS w eff = for w d = 

J J x-\-y+z 

-0.7,6 = 0.5. 


following form: 

—— = — 3x(w d — 2(w d x + w r z) j — 96a;z, (30) 

dy 

= 6y{w d x + w r z) + 9bz(x-y), 

= —3z(w r — 2(w d x + w r z)^+9byz, 

Critical points and the eigen values are obtained as the 
following: 
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• A 4 : {i = 0 ,y = y,z = 0}, A, = {0,-1, -3w d }. 
This is an unstable point. 

• B 4 : {x = 0,y = 0,z = 1/2}, A/ = {1,1 — f, 1 — 
// — 3 w d }. This point is also unstable. 

• C 4 : {x = l/2,y = 0 ,z = 0},Aj = {3w d - 
1, 3 Wd, 3 w d }. The point here is stable. 


• £4: 


J_ _ _ _ Wd 

~ 96’ V ~ 96’ " ~ 36 


(31) 


we have a cubic equation for determining eigen¬ 
values as a function of the absolute cosmological 
parameters { 6 , w d }: 


A 3 _ 3wjX2 _ 2w d [-l + 3w,) x 
60 

| w d 2 (—4 + 96 + 6 w d ) _ 0 


(32) 


This cubic Eq. can be solved by quadratures. For 
6 > 0 , Wd < 0 , it is possible to numerically show 
that one eigenvalue always is positive A 3 > 0. So 
the system is unstable. 



FIG. 10: Model IV: Solutions x(N),y{N), z(N) for w d = 
—0.7,6 = 0.5. Solid x(N), dot y(N) and dash z(N) 


• £4: 


1 2 — 96 — Qw d 
X = “18 6 ’ 

27w d b + 18w^ — 96 — 12 w d + 2 
V = -186 + 816 2 
^ _ -6w% + 2 w d 

Z ~ —66 + 2 76 2 ' 

eigen values: 


(33) 


A d - 


18w d (—1 + 36 + w d ) 
-2 + 96 


X 2 +AX + B = 0. (34) 


here 

3 Aw- 1 (-2 + 9 6) 2 6 = -8 - 324 6 2 - 72 w d 2 + 72 6 (35) 

+48 w d — 972 w d b 2 y + 4374 w d b 3 y + 972 w d 3 b + 1458 w d 2 b 2 
-648 w d 2 b - 486 w d b 2 - 1458 b 3 y + 324 b 2 y + 729 6 3 



-B (-2 + 9 6) 2 6 (-1 + 3 Wd)- 1 w d ~ 2 
= —108 6 + 8 — 216 w d by + 972 w d b 2 y + 216 w d 3 
+324 w d 2 b — 486 w d b 2 


(36) 


FIG. 11: Model IV: Phase space for Wd 
shows an attractor behavior. 


-0.7,6 = 0.5. It 


+1458 b 3 y - 144 w d 2 + 108 w d b - 648 6 2 y 

+72 by - 729 6 3 + 48 6 6 2 L CONCLUSION 


The eigenvalues has at least one positive real num¬ 
ber A + > 0. So the system is still unstable. 


In Fig. (10) a numerical solution is developed for density 
function(s). All types of the densities are decaying. 

An attractor which was started from an initial point is 
observed in Fig. ©■ 

An affective EoS with range w e ff < —0.5 is drawn in 
Fig. ©. 


The final achievement of modified gravity in this di¬ 
rection was the extension of the method of the disformal 
transformations of the metric of arbitrary dimensions, 
successfully used by authors in the investigation of cos¬ 
mological as well as of gravitational equalities, to any 
system whatever of mutually interacting components like 
dark energy and dark matter. We further supposed that 
the dark components are interacting like fluids, but each 
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N 


FIG. 12: Model IV: Effective EoS w eff = for w d = 

—0.7, b = 0.5. 


perceiving what was passing in the other, and acting in 
consequence by interaction term I\. This letter was the 
first to investigate and describe in mimetic gravity the 
fact that a cosmological solution at a distance had the 
power of invariance from initial condition to final state, 
and it also found that in some interacting cases the stable 
submanifold was produced. 
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